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Recent controversy on the quantum dots dephasing mechanisms (between pure and inelastic) is 
re-examined by isolating the quantum dots from their substrate by using the appropriate limits of 
the ionization energy theory and the quantum adiabatic theorem. When the phonons in the quan- 
tum dots are isolated adiabatically from the phonons in the substrate, the elastic or pure dephasing 
becomes the dominant mechanism. On the other hand, for the case where the phonons from the sub- 
strate are non-adiabatically coupled to the quantum dots, the inelastic dephasing process takes over. 
This switch-over is due to different elemental composition in quantum dots as compared to its sub- 
strate. We also provide unambiguous analyses as to understand why GaAs/AlGaAs quantum dots 
may only have pure dephasing while InAs/GaAs quantum dots give rise to the inelastic dephasing 
as the dominant mechanism. Our study accentuates the importance of the elemental composition 
(of both quantum dots and substrate) in evaluating the dephasing mechanisms of quantum dots. 
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I. INTRODUCTION 

The electron-phonon (e-ph) interaction in quantum 
dots (QDs) is evaluated theoretically by invoking the 
discontinuous dielectric property (between QDs and sub- 
strate) and quantum adiabatic approximation in order to 
understand the coupling between electrons in QDs and 
the phonons in both QDs and substrate. In order to do 
so, we propose here that the elemental composition of the 
QDs and the substrate need to be explicitly considered. 
We further show that such consideration will lead us to 
understand the dephasing mechanisms in QDs accurately 
for any non-frcc-electron QDs and substrate. Apparently, 
(1) the e-ph interaction within the QDs and (2) between 
QDs and the substrate are crucial to the understanding 
of the dephasing mechanisms Pure dephasing pro- 
cess is defined as due to the elastic interaction between 
electrons and phonons that only further corrects the un- 
perturbed energy levels in QDs. Such a process does not 
excite the electrons or change the excited carrier density 
or the carrier occupation numbers 0, H[. Whereas, in- 
clastic dephasing requires phonons to assist the electrons 
and/or holes to be relaxed to a lower energy-level, or elec- 
tron relaxation due to phonon emission [I, i, Q. Both 
pure- and inelastic-dephasing processes hinder the stor- 
age of quantum information, which are characterized by 
the timescales, Ti and T2, respectively 0]. The respec- 
tive T\ and T2 timescales are also known as the phase 
and population relaxation lifetimes. 

Presently, there are several theories on these dephasing 
processes and are given in the Refs. 0, d, 0, [H| that 
employ the adiabatic approximation. In those treatments 
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however, the adiabaticity is invoked between the nuclear 
motion and the electronic excitation with respect to de- 
generacies .without incorporating the phonons from the 
substrate 0, [1, d, QjJ LIU- Apart from that, a polaronic 
model based on the polaron relaxation was proposed to 
understand the relaxation mechanisms in QDs [12j , which 
employs the app roach of Klcmens, Barman and Srivas- 
tava [ID, [ID, Gil- In their study, only Frohlich polarons 
or the electron-LO (longitudinal optical) phonon interac- 
tion was considered. This model is an improved ver- 
sion of other reported relaxation models, namely, (1) the 
semi-classical approach [l6[ and (2) a phenomenological 
model that invoked the anharmonic mechanism for the 
bulk LO phonons within the Fermi golden rule frame- 
work [It], [HI, HH. In all these approaches, (1) the origin 
of phonons is from the QDs, where the LO phonon contri- 
butions from the substrate or wetting layer are ignored 
by assuming the substrate phonons are the low-energy 
acoustic phonons, and (2) the phonons are independent 
of elemental composition. 

Here, we take another step forward to understand the 
influence of substrate phonons (both acoustic and opti- 
cal) on QDs dephasing mechanisms (both pure and in- 
elastic) in detail. The influence of substrate phonons will 
be captured by studying the elemental compositions in 
QD and substrate materials separately. We will invoke 
the ionization energy theory (IET) to analyze how dif- 
ferent elemental compositions in QDs and substrate will 
give rise to pure or inelastic dephasing processes or both. 
The reasons for incorporating the elemental compositions 
in QDs and substrate are (1) to investigate the reports 
that claim the matrix and/or substrate are indeed influ- 
encin g bo th the electronic and phononic properties of the 
QDs (2O, [HI, IH, HI Hl|, and (2) to develop a strategy 
based on the IET to distinguish the contradicting results 
obtained by Sanguinetti et al. Q and others: Zibik et 
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al. 0, [H| and Chernikov et al. [26| . For example, pure 
dephasing was found to be responsible in GaAs QDs on 
Alo.3Gao.7As substrate Q. Whereas, inelastic dephasing 
played the major role in InAs QDs embedded in the GaAs 
matrix [El I25I [26| . Therefore, our primary intention here 
is to understand why and how different constituent atoms 
in QDs and substrate materials may contribute to the dif- 
ferent electronic relaxation or dephasing mechanisms in 
QDs. 

In order to achieve this, we will employ the IET, which 
has been developed earlier [U [H, [2^, [3(| ■ The reason 
to use the IET is because the theory is straightforward 
and the physical mechanisms derived from the IET can 
be directly related to the atomic constituents of any non- 
free-electron systems [3l|, HI, H3| • We organize the paper 
in the following order. A brief introduction to the IET 
is given in the following section. Subsequently, techni- 
cal discussion are developed on the technique of isolating 
the QDs from its substrate. The QDs are first electron- 
ically isolated and followed by adiabatic isolation. De- 
tailed discussion on experimental proofs and the possible 
applications of the IET arc highlighted with predictions. 
In addition, wc also explain how the different dephas- 
ing rates due to defects, impurities, electron-electron and 
electron-phonon interactions, and spins are related to the 
ionization energy concept. 



II. ELEMENTS OF THE IONIZATION ENERGY 
THEORY 

We start from the many-body Hamiltonian, which is 
given by [H HI 



H<p=(Eo±Z)<p, 



(1) 



where, the eigenvalue is exactly equals to Eo±£. Here, 
E is the total energy of the system at zero temperature 
(T = K), £ is the energy- level difference in a given 
atom or the ionization energy. Importantly, the IET pre- 
sented here can be used to study any non-free-electron 
materials. If a given material is a free-electron metal, 
then Eq. ([1]) reduces to the standard time-independent 
Schrodinger equation given by, Hip — Etp, where one 
needs other theoretical and computational methods to 
solve it by means of variational principle [34j . For solids 
however, the parameter £ in the total energy, Eq±^ refers 
to the energy level differences in solids, which is difficult 
to be determined. Therefore, wc will use the ionization 
energy approximation, £ oc £7, where Ej is the atomic 
ionization energy (for an isolated or free atom). Math- 
ematically, the ionization energy approximation can be 
written as 



E ±^E ±Y,— > 



where, the subscript i counts the first, second, z 
ionization energy of each constituent atom for a given 
material. Here, J^i En/ z gives the changes to the aver- 
age ionization energy of a given system. In addition, in 
this approximation, one needs to rely on the accuracy of 
the ionic valence states of the constituent atoms and as 
such, knowing accurate valence state values are impor- 
tant, and usually, they are easily predictable. For exam- 
ple, these valence states can be obtained from their stable 
oxidation states and these states are also known to vary 
significantly as a result of defect-formation and/or other 
structural deformations Subsequently, we can sub- 
stitute the new total energy, Eq ± Ei (after the approx- 
imation) into the ionization energy based Fermi-Dirac 
statistics [27], [29| as given below 



f e (Eo,Ei) 



e [{E +Ei)-E^>]/k B T + 1 

fh(Eo,Ex) = , 0) - * . (3) 



-{E -E,)]/k B T + 1 



Note here that for QDs, E^ can be regarded as the 
highest occupied energy-level at T = K, and fcg is 
the Boltzmann constant. Furthermore, f e (Eo,Ei) and 
fh(E ,Ei) are the probability functions for electrons and 
holes, respectively Even though Eq. ([T]) is for T = 
K, the temperature effect can be taken into account via 
Eq. ([3]) . Note here that the appearance of £ in any equa- 
tions simply means that the ionization energy approxi- 
mation has not been invoked. After applying the approx- 
imation, one must replace the £ with Ej in all the equa- 
tions. This approximation will be used later to calculate 
the average ionization energy of the elements present in 
a given QD. The size effect is easily obtained from a sim- 
ple one-dimensional infinite square-well potential, which 
is given by [HI 



y/2m(E ± Ca)r 



(4) 



(2) 



Here, a is the width of the potential well and wc can 
see the inverse proportionality between a and £, while 
to, h and n denote the electronic mass, Planck constant 
and the principal quantum number, respectively. Impor- 
tantly, Ei is unique for each atom that will assist us to 
capture the effect of different constituent atoms in QD 
and substrates on their electronic excitation probability 
and the electron-phonon interaction. Therefore, the only 
input parameter is the type of atoms or ions that may 
exist in a given sample, and their appropriate valence 
states. More details about the IET and its approxima- 
tion can be found in Refs. [2?], [H, [2iJ [3(j . In the subse- 
quent sections however, we will use £ throughout instead 
of Ei for notational clarity. Thus, any values calculated 
for £ implies that we have invoked the ionization energy 
approximation. 
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III. ELECTRONIC AND ADIABATIC 
ISOLATION OF QDS FROM THE SUBSTRATE 

In the following sections, the electronic interaction be- 
tween the arrays of QDs and the substrate is evaluated 
with respect to the possibility of isolating the electronic 
property of QDs from the substrate. After that, these 
QDs arc separately (1) coupled to QD and substrate 
phonons, and (2) isolated from those phonons adiabat- 
ically. This adiabatic approach is used to understand 
the influence of the substrate and QD phonons with QD 
electrons. We also compare our theoretical results with 
recent experimental observations. 

A. Electronically isolated QDs from the substrate 

Let us first evaluate the charge distribution of the QDs 
with respect to its substrate using the Gauss and Green's 
theorems. From the Gauss theorem [3a |. we can write 



writing Eq. (J6]) is to obtain separate charge distributions 
in QDs and substrate. By isolating them, one can further 
evaluate the strength of charge accumulation in QDs, as 
compared to their substrate. As such, one can rewrite 
Eq. ([6]) to arrive at 

= [ [v ■ 1 SUB - v • A^ B ]dn 

Jn 

= [ [V ■ (e SUB W SUB ) - V • ( f « D V^ D )]<l!l 
Jn 

= [ [ £ SUB V 2^SUB _ e QD V 2 y QD ]dQ; (?) 

Jn 

where, dSj is an element of the substrate surface, and 
dS p is for QD surface. For QD arrays on a substrate, we 
will consider three conditions, namely, 



E 



A ■ tij dS 1 ; = / V • A dCl, 



m surface 



Substrate 



q-1 surface 



m-1 surface 



q^ 1 surface 



(5) 



FIG. 1: Two-dimensional schematic representation of an iso- 
lated QD and substrate surfaces with m — 1 and q — 1 closed 
surfaces, enclosed by the m th and q i surfaces, respectively. 



t SUB ^ ^QD 
t SUB 



condition (i), 
condition (ii), 
condition (iii). 



These conditions will be used to identify whether the 
electrons in QDs are both electronically and/or adiabat- 
ically isolated from their substrate. Satisfying condition 
(i) implies that one can use the last integral in Eq. ([7]) to 
obtain (after using Poisson equation, eV 2 U = —p/eo) 



/ \SUB 
H\ induced dn 



< 


/ 




Jn 



P( r )induccd 



dQ 



eo 



(8) 



where, f2, rij and A denote the volume, unit vector 
normal to the surface and any given vector, respectively 
with j = 1, 2,... m counts the number of closed surfaces 
(see Fig. [T]). We can now use Eq. (O by first defining, 
A = D = eE = -e\7V, V, D and E denote the elec- 
trostatic potential, charge displacement and the electric 
field, respectively. Subsequently, by substituting them 
into Eq. ^ we can obtain 



E 



eW 



rijdSj 



J 

Jn 



V • eVVdn 



eV 2 Vdn. (6) 



Thus far, we have not explicitly described the quantity, 
e, which is a dielectric constant. However, e is obviously 
discontinuous at the boundary between the substrate and 
QDs. For a given array of QDs on a particular substrate, 
there will be m — 1 surfaces enclosed by the m surface 
and similarly, there will be q — 1 surfaces enclosed by 
the 9 th surface for the substrate (SUB). The purpose of 



where, p(r) is the charge density, which is given by [2t 



P(r) 



Po(r) +p(r) ind 
3noe 2 



= n e + 



2E° F 



V(r)e 



(9) 



where e is the electronic charge, while A = 
(127reo/e 2 )as, in which is the Bohr radius and eo 
denotes the permittivity of free space. Note here that 
Po(t) and hq are the charge and carrier densities, re- 
spectively at T — K and without any external dis- 
turbances. Whereas, the induced charge density, p(r)i n a 
capture the charge accumulation due to temperature, 
T > K and due to other external disturbances. The 
potential, V(r) = V(r) cxt + U(r) ind , where U(r) cxt is the 
external potential. If p(r)i n( j = 0, then we will arrive at 
the Green's reciprocation theorem [3{| . For example, af- 
ter letting p = j and q = m, so that, f2 SUB = fJ QD (for 
mathematical convenience) 
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E 



3=1 



3=1 J ^ 

or equivalently 



dSi « 0, 



(10) 



E 

3 = 1 



,/.S, - / 7j QD dS, 

5 7 " Sj 



0, (11) 



where 7 denotes the surface charge density and subse- 
quently, we can arrive at 



E e ?° 

3=1 



E e f B > 
3=1 



because, 7 



SUB 



7 



QD 



or Po UB ( r ) 



(12) 



p^ D (r), which is to 



say, both QDs and the substrate are approximately un- 
charged in the absence of temperature and external dis- 
turbances. Now, we can repeat this procedure exactly by 
only switching the indices, QD and SUB to obtain condi- 
tion (ii), £ SUB < £ QD . Hence, both conditions (i) and (ii) 
imply that the QDs are electronically isolated from the 
substrate. However, condition (hi), £Q D w £ SUB gives 
rise to e-e interaction between QDs and its substrate. 
For example, from condition (i), we have <C £ SUB — > 
w qd <<; w sub _> f^Ejfiv » f e (E ,Ej) SVB , thus, 
the characteristic electronic (subscript e) excitation prob- 
ability for QDs is much higher compared to its substrate, 
and vice versa for condition (ii) (for example, see Eq. ([3])). 



SUB 



Note here that ^ D /h = uj^ d and £ SUB //i ; 

Where condition (iii) holds, the electronic properties 
of the substrate will also affect QDs and vice versa. This 
means that the electronic excitation probability for both 
QDs and the substrate are identical and therefore they 
(QDs and the substrate) are electronically coupled. In 
other words, condition (iii) also implies that the confine- 
ment energy is close to zero. However, this electronic 
coupling does not occur in real QDs that have finite con- 
finement energies. This may explain why the QDs are 
always electronically isolated if the confinement energy ^ 
0. Shapes of QD structures also significantly affect their 
electronic properties [HI, H3, [H, [33] • This shape depen- 
dence will give rise to anisotropic ionization energy, and 
becomes important, if and only if we are comparing QDs 
of different shapes. In reality, for a given nanoscale syn- 
thesis process, shapes of QDs are quite similar, whereas 
size non- uniformities are quite substantial [13, Sl|. In 
cases where QDs have different shapes, £ can be spa- 
tially averaged by incorporating the QD size, a correctly 
(will be a complicated algebraic function) in Eq. 



B. Electron-phonon interaction and 
adiabatic-isolation of QDs from the substrate 

The next issue here is to understand the phonon contri- 
bution to the electronic properties of QDs. We will inves- 
tigate the possibility of isolating the substrate phonons 
from QDs, or coupling them to QDs. The motivation to 
study these electronic (given earlier) and adiabatic iso- 
lations are important so as to evaluate the dephasing 
mechanisms namely, elastic (pure) and inelastic dephas- 
ing. We will address both of these mechanisms separately 
here with the IET and the quantum adiabatic approxi- 
mation. It has been shown in our earlier work, that the 
phonon frequencies of both optical and acoustic, for ID 
diatomic system with two different ions (masses, Mi and 
M 2 ) is given by (27j] 



e A(£-E°) 



[A±B*], 



(13) 



2M X M 2 
A = (Q + G)(M 1 +M 2 ), 
B = {Q + G) 2 (M 1 +M 2 ) 2 - 
4MiM 2 [(Q + Gf -Q 2 -G 2 - 2QG cos(ku)] . 

where, Q and G are the interaction potential constants, 
ii is a lattice site, whereas, the screened e-ph Coulomb 
potential is given by [13] 



Uep(k,k*) 



1 

fie 



|k-k*|2+i^exp[A(£° -£)] 



,(14) 



which is in exact form, with the screened e-e Coulomb 
potential [29| . 



Uee(k) 



1 

£le 



k 2 +x2 C xp[A(^-0] 



(15) 



where, k — k* = q is due to the crystal momentum 
conservation. Here, k and K s denote the wavevector and 
Thomas-Fermi wavenumber, respectively. Figure [5] shows 
the V cp dependence on the ionization energy and the QD 



size, a that follows from Eqs. (Ti"4"|) and 0, while the 
inset captures the definition of the ionization energy in 
a given QD. Note here that the ionization energy here 
implies the energy-level difference between the ground 
and excited states. 

From Eq. (fl3|) . condition (i) implies wjj <C Wp™ — > 
t QD ^sub^ w h er6j t QD an( j ^sub denote the character- 
istic timescales (not the relaxation lifetimes or dephasing 
rates) for the phonons in QDs and substrate, respectively. 
Note here that w p h can be due to optical or acoustic or 
both. 



A If Q D ^ SUB 

Analysis jor w^ h >C w ph 



Claim 1: The electrons in QDs can be made adiabat- 
ically coupled to the substrate or QD phonons or both 
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5 10 15 20 

QD size, a (arb. units) 



FIG. 2: Electron-phonon interaction energy as a function 
of QD size, a for two different ionization energy values, 
Eil < Ei2. These plots were obtained from Eqs. ([T4)l and Q. 
The INSET shows the definition of the ionization energy in a 
QD, with appropriate labels for the ground and excited state 
energies. Note here that the ionization energy is the energy- 
level difference between these two states. 



by Zibik ct al. [H[. 

For example, the inequality, C3 (and/or D4) have been 
achieved recently by deliberately changing the elemental 
composition of InGaAs/GaAs QDs via annealing [25] ]. 
Note here that ^° aAs < ^As- They have (a) ob- 
taincd the inelastic dephasing due to C3 and/or D4 sat- 
isfying the theoretical results discussed above, and (b) 
achieved longer relaxation lifetime in QDs (r^ D ) by en- 
forcing £Q D < fi,K7p™ ,( ^ D . In other words, they system- 
atically reduce the photoluminescence (PL) energy peaks 
by annealing their samples. Reducing PL energy peaks 
means decreasing £Q D and therefore increasing the mag- 
nitude of fiWp® — £Q D . As a result, they have enforced the 
inequalities, < r^ D and/or t SVB < t^ d as predicted 
by the IET and condition (i). 



2. Analysis for £ QD > £ : 



by employing condition (i). 

Proof 1: From condition (i), we have £ QD < £ SUB 
that implies « w 5 ™ because £ QD ' SUB oc w^ SVB 
(from Eq. (fT5| ). The next step is to identify the 
relationship between and Wp h D,SUB , where we can 
define four inequalities, 

Al ^ D >^SUB^ e QD_^SUB^^UB >r QD ; 

B2 e QD > e QD -p/i QD =► i QD > rQ°, 

C3 C QD <^ph B ^ eQD -P /lSUB ^ iSUB <^ QD I 
D4 e QD <^ h D ^eQ D -p/ l Q D ^tQ D <rQ D , 

where r^ D is the electronic relaxation lifetime or de- 
phasing rate in the QDs. Both Al and B2 will give rise to 
pure dephasing because the phonons timescales from the 
substrate (i SUB ) and QDs (^ D ) are low energy phonons 
compared to the QD energy-level difference (£^ D ). In- 
voking condition (i) implies Al and B2 are unlikely and 
thus pure dephasing is negligible. On the other hand, the 
inequalities in C3 and D4 satisfy condition (i) and will 
lead us to inelastic dephasing because of the adiabatic 
coupling between the electrons in the QDs and phonons 
from both QDs and substrate. Hence, there is always a 
finite probability for the absorption and emission of QD 
(e QD -p/i QD ) and/or substrate (e QD -p/i SUB ) phonons ■ 

Now, we know that condition (i) gives rise to inelastic 
dephasing as the dominant mechanism. Therefore, the 
relaxation lifetime, r^ D can only be made longer by fur- 
ther enforcing the inequalities < r^ D and/or i SUB < 
T QD j n th er words, we need to satisfy and enforce the 
inequalities, £ QD < hw^ and/or £ QD < hw^ B . For 
example, enforcing the inequality, ^ D < tvw^ means 
fvuS^ — £Q D — > oo. This is exactly what have been done 



Claim 2: The electrons in QDs can also be made adia- 
batically independent from the substrate phonons by em- 
ploying condition (ii), which is ^ D ^> £ SUB . In this case, 
we have w^ h D > w^ B -» i QD < i SUB and therefore, 
the e^ D -p/i SUB interaction between QD electrons and the 
substrate phonons can be assumed to be switched-off (no 
inelastic dephasing). To see this conclusion clearly, we 
need to recall proof 1 given earlier. 

Proof 2: After invoking condition (ii), C3 and D4 are 
unlikely and therefore the inelastic dephasing is negligi- 
ble. Consequently, we are left with the inequalities given 
in Al and B2, which satisfy condition (ii) and pure de- 
phasing as the dominant mechanism ■ 

Here, the relaxation lifetime, r^ can only be made 
longer by further enforcing the inequalities i SUB > r^ 
and/or > r^ D as given in Al and B2. Based on 
condition (ii), we need to enforce the inequalities, £Q D > 
hw^ and/or ^ D > hw^ B with appropriate changes 
to the elemental composition of the QDs or substrate or 
both. 



3. Analysis for » B 

Finally, if the QDs were to satisfy condition (iii) , then 
these QDs are just dots with no electronic confinement 
effect and however, their e-ph interaction cannot be un- 
derstood without additional complications. In this case, 
the ionic mass will start to play a crucial role as given in 
Eq. (|13p . In other words, condition (iii) does not guar- 
antee the relation, = w^ B , and the reduced mass 
of the QDs and substrate will determine the correct rela- 
tionship. For example, the reduced mass (M r ) is defined 
as 1/M r = J^jl/Mj and even if condition (iii) is ful- 
filled, we may still have Af r QD ^ M^ UB , which will lead 
t° w ph ^ ^ph 8 - To understand this effect clearly, let 
us first assume that the QDs and the substrate are of 
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same materials and elemental composition, and we do 
not intend to change the elemental composition. This 
assumption implies, M^ D = A/^ UB and the inequality, 
t'p h D > w^ B is due to £ QD > £ SUB , which is entirely due 



to size. In addition, in this hypothetical case, condition 
(iii) can be achieved by making the size of the of QDs 
larger. By doing so, these QDs lose both their electronic 
and phononic confinements. 

On the other hand, if the QDs and the substrate ma- 
terials are different (M r QD ^ M r SUB ), then condition (iii) 
can be achieved by changing the elemental composition of 
the QDs. As a consequence, we have both M r QD ^ M^ UB 
and £ QD = £ SUB , which imply that the QDs are not elec- 
tronically isolated, but are still isolated with respect to 
phonons. In this second scenario, it is essential to note 
that one cannot vary £ QD independent of AI^ U . Apart 
from that, small changes in the reduced mass does not 
play a significant role in the presence of condition (i) and 
(ii) because the strength of the electronic polarizability 
for each ion is more effective in determining the w p h- This 
electronic polarizability is related to the interaction po- 
tential constants, Qe x ^~ E ^ and Ge x ^~ Ep \ which takes 
the deformable-ion effect into account. Therefore, if we 
have a situation where M r QD > M r SUB -> u^ h D < w^ B 

and C QD > £ SUB -> ttfS? > w^ B , then M r and £ com- 



pete with each other. In this case, £ is more sensitive than 
M r because the ionization energy that control the rigidity 
(or deformability) of the ions determines the u> p h- For ex- 
ample, slight changes in the reduced mass can change the 
phonon frequency, but not as effective as the ionization 
energy that control the ions deformability. However, any 
huge changes to M r will significantly reduce the effect of 
the ionization energy on u> p h- 



IV. EXPERIMENTAL PROOFS AND 
APPLICATIONS 



In this section we will scrutinize QDs made of two dif- 
ferent materials in which, one of them has been studied 
by Sanguinetti et al. |2j while the other one was investi- 
gated by Zibik et al. 0] . Note that the QDs and their re- 
spective substrates discussed here have Zincblende crys- 
tal structure [H| and therefore, Eqs. (TT3]), (0 and $To§ 
can be applied directly [27| . 



A. GaAs QDs on Alo.3Gao.7As substrate 



QD sample, is much larger than the calculated activation 
energy (« 30 meV, related to GaAs longitudinal optical 
(LO) phonon). Therefore, the LO phonons do not con- 
tribute to inelastic dephasing As such, pure dephas- 
ing is the likely process with GaAs LO phonons. These 
low-energy LO phonons do not cause electron relaxation 
from an excited state to the ground state. 

Let us now apply our theory described earlier using 
conditions (i) and (ii) to GaAs QDs on Alo.3Gao.7As sub- 
strate. The average ionization energies for both Ga and 
As are given by, £ Ga 3+ = (578.8 + 1979.3 + 2963)/3 = 
1840 kJmol" 1 and £ As3+ = (947 + 1798 + 2735)/3 = 
1827 kJmol -1 , respectively. Prior to averaging, all the 
ionization energies mentioned in this work were obtained 
from Ref. 0]. This gives the total value of ^g T = 1827 
+ 1840 = 3667 kJmol" 1 . Similarly, the total average for 
the substrate can be calculated as, £ca 3 + = [0.7 x (578.8 + 
1979.3 + 2963)]/3 = 1288 kJmol" 1 , £ A1 3+ = [0.3x(577.5 
+ 1816.7 + 2744.8)]/3 = 514 kJmol" 1 , therefore, the to- 
tal value for the substrate is = 1288 + 1827 + 514 
= 3629 kJmol" 1 . Thus, the difference is, £™t - £™t 

= 3667 - 3629 = 38 kJmor 1 = 400 meV atom" 1 and 
we have £or < ^tot ^ w sub < W QD _ ^ub > 4 qd 

Obviously, this satisfies condition (ii) and gives rise to 
adiabatically independent QDs and elastic dephasing is 
expected from our theory, in accordance with the exper- 
imental results and interpretations given in Ref. Q. 

Apart from that, the inequality, £,s\ib < £qd T implies 
that the electrons in the substrate has a higher proba- 
bility to be excited at any given temperature, compared 
to the electrons in QDs. As a consequence, the substrate 
ions deformability is larger due to large electronic excita- 
tion probability, which in turn, explains why the lattice 
vibrational frequency for the QDs is larger than the sub- 
strate [I?}. This means that, the ions in the QDs arc 
more rigid compared to their substrate. 

Now, using the reduced mass (in atomic mass unit) as 
the main argument, we can obtain M^ UB (28.98) < M r QD 
(36.11) -> w^ B > wjf, which is in contradiction with 
the analysis stated above. However, as explained earlier, 
we have assumed that the effect of the small reduced 
mass is not as substantial as the ionization energy effect 
in accordance with Ref. [27j • Physically, what this means 
is that when the change in the reduced mass competes 
with the change in the ionization energy, then for as long 
as the change in the reduced mass is relatively small, the 
ionization energy effect will stay significant. 



The former group [2j studied GaAs QDs on 
Alo.3Gao.7As substrate grown using modified-droplet epi- 
taxy with QD density of the order of 10 8 cm -2 . They 
observed that the exciton dephasing of the QDs is in- 
dependent of confinement energy with no quantum-size 
effect. This conclusion was based on their temperature- 
dependent PL spectra measurements of which, the 
confinement-energy (as 160 meV) of one of their smaller 



B. InAs QDs in GaAs matrix 

The second evidence comes from the work of Zibik et 
al. 0. Their QDs are made from InAs grown layer- 
by-layer with GaAs as the matrix and barrier. They 
used the four- wave-mixing (FWM) spectroscopy to study 
the dephasing mechanism in QDs. The samples were 
grown using the molecular beam epitaxy in the Stranski- 
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Krashtanow mode and the QDs were separated by 50 nm 
wide GaAs barriers to prevent structural and electronic 
coupling between the QD layers [f|. Both their calcula- 
tions and experimental results point toward intersublevel 
(inelastic) dephasing in electron-doped QDs, which re- 
vealed oscillatory behavior for the polarization decay (for 
times < 5 ps). By repeating our £ averaging as shown 
above, and after identifying the GaAs matrix/barrier as 
the substrate, we can arrive at ^svb ~ £qd T = 3667 
- 3521 = 146 kJmor 1 = 1500 meV atom" 1 . In this 



< 



case, we have > Cqd 1 

£Q D . Hence, InAs QDs satisfy condition (i) in which, the 
phonons from GaAs coupled to QDs non-adiabatically 
and favor the inelastic dephasing mechanism. 

The inequality, £guB > £q§ T a l so implies that the 
electrons in the substrate have a lower probability to 
be excited at any given temperature, compared to the 
electrons in QDs. As a consequence, the substrate ions 
deformability is smaller due to smaller electronic excita- 
tion probability, which in turn, explains why the lattice 
vibrational frequency for the QDs is smaller than the sub- 
strate [27|. This means that, the ions in the QDs are less 
rigid compared to their substrate. As anticipated, our re- 
sults again agree qualitatively with the results calculated 
and measured experimentally in Refs. p|. 

Unlike the GaAs QDs discussed earlier, the reduced 
mass effect for InAs QDs also further supports the conclu- 
sion that > tift . In other words, for InAs/GaAs 



system, we find that M r SUB (36.11) < M r QD (45.34) 



„SUB 



> W 



QD 
ph 



which eventually satisfies the inequality, 

t SUB < t^ u . One can surmise here that there is no com- 
petition between the ionization energy effect and the re- 
duced mass contribution to the phonons. 



C. Further analysis 



system, which can be obtained from (using Eq. (8) given 
in Ref. (U) 



Eq ± £ = -Ekinetic + Vcoulomb + V b ™dy Y 

L J y 



(16) 



where 



/-rrmanyv 
Nobody / 

Ei ' 



(17) 



Consequently, we can identify (V b ™^ y y ) as the atomic 
screened Coulomb potential by assuming that QDs are 
artificial atoms, which is given by (2?| (we apply Eq. (TT|) 
to two-electron hydrogen-like atomic systems) 



H = H + V SC , 
and its solution is given by [13] 



(18) 



(if) = 2Z Ei + (V) ac = 2Z E\ + j— 1 —. (19) 



[2Z + a B (r} 5 ' 



Where, Z is the atomic number, 



y _ e e -nre\ x{ 5) _ e p-QT 



47reor 



(20) 



Here, /i is the screening parameter's constant of pro- 
portionality and 



At this juncture, one may wonder two possibilities, 
namely, (1) if the size of the QDs arc much smaller than 
the substrate, then the inequality, > £qd T as stated 

for InAs QDs may not be true for all QD sizes, and as 
a result of this, (2) there is a possibility to grow dots 
that are not electronically isolated, in other words, during 
growth, the dots at a certain small size may satisfy Csub" 
= £q§ T , or even tg$g < £^g T for very small InAs QDs. 
Hence, one can indeed in principle grow dots with zero 
confinement energy by choosing the appropriate QD size 
and elemental composition for a given substrate. How- 
ever, the changes in ionization energy due to size is much 
less compared to the change due to elemental compo- 
sition (between InAs QDs and GaAs matrix). In other 
words, the inequality, £gfj^ > £qd T will not be physically 
reversible for InAs QDs surrounded by GaAs matrix by 
only changing the QD size because of large contribution 
due to elemental composition (1.5 eV atom -1 ). To see 
this effect clearly, let us write down the respective ioniza- 
tion energy approximation equations for the InAs/GaAs 



E x 



2h 2 \4ire 



(21) 



Using Eqs. (TT5|) and (TiU)) , we can surmise that indeed 
smaller £tmik gives rise to smaller confinement. Now, the 
change in £buik due to elemental composition is ^uik = 
Cbuik — £bu?k- Whereas, the change in the £ due to size 
is ^Le = ££f e B ~~ €s&e- Since the change in the ionization 
energy due to size is smaller compared to elemental com- 
position, we can assume that l^uik > l^sizel' which was 
done earlier. Even if |^ ulk | > \8^ izc \ is not true, then it is 
possible to grow very small QDs that are electronically 
coupled to the substrate. Recall here that this electronic 
coupling means that the electronic excitation probabili- 
ties for both QDs and the substrate are identical. 

In summary, using the IET, it is possible for us to select 
suitable substrate and QD materials (elemental composi- 
tion) in order to switch-off the contribution of substrate 
phonons. This will help the experimenters to grow QDs 
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FIG. 3: The evolution of the QD's ionization en- 
ergy with respect to systematic substitutional doping for 
Ga2-a;Asa;/Alo.3Gao.7As and h^-^As^/GaAs (Inset). Sub- 
stituting atomic Ga with As decreases the total ionization 
energy of the QDs (Ga2- a: As : i;) but it is still larger than the 
total ionization energy of the substrate (Alo.3Gao.7As), which 
is 3629 kJmoP 1 (follow the arrow). INSET: The total ion- 
ization energy for the Ina-aAsa; QDs increases with x but it 
will still be smaller than the total ionization energy of the 
substrate (3667 kJmol -1 , follow the arrow), regardless of x. 
Therefore, both QDs are electronically isolated from their re- 
spective substrates for all doping, x. See text for details. 



with specific application in mind, be it for quantum com- 
puting (with negligible phonon contribution from sub- 
strate) or photovoltaic cells (with large confinement en- 
ergy and small e-ph scattering). 

For example, Fig. [3] indicates what to expect if we were 
to change the doping in G&2- x As x and h^-^As^ QDs 
systematically. Regardless of doping x, we will find that 
the total ionization energy of the Ga2- x As x QDs will 
always be larger than their substrate's value. The rea- 
son is that the total ionization energies for both x = 2 
and x = are still larger than their substrate's averaged 
value, 3629 kJmol -1 . On the other hand, the total ioniza- 
tion energies for In2_ x As x QDs for both x = 2 and x = 
are smaller than their substrate's value, which is 3667 
kJmol -1 (see the inset of Fig. [3j> - Consequently, these 
QDs will always be electronically isolated regardless of 
the doping x. In other words, Ga2_ x As x QDs will re- 
main adiabatically independent from their substrate for 
all doping x (pure dephasing). Whereas, for In2_ x As x 
QDs, the substrate phonons will keep influencing the elec- 
trons in the QDs for all doping x (inelastic dephasing). 
It is important to note here that even if we were to in- 
crease the size of these QDs to form thin films, these thin 
films will still be electronically isolated from their sub- 
strate due to their different elemental composition. That 
is why thin films of different elemental composition from 
their substrate will have different electronic properties 
compared to their substrate. This will remain true for as 
long as the elemental composition between the substrate 
and the thin films are different. 



V. DEPHASING MECHANISMS AND 
IONIZATION ENERGY 

Here, we discuss the various types of dephasing rates 
(relaxation lifetimes, r) due to defects (def), impurities 
(imp: magnetic and non-magnetic), spin, e-e and e-ph 
interactions, and how they are related to the IET. These 
effects arc first analyzed in the absence of any external 
perturbations such as laser (to create excitons) and ap- 
plied magnetic fields. Subsequently, further evaluations 
are carried out based on the changes in the spin transition 
probability by identifying the fact that the inelastic de- 
phasing is favorable with increasing transition probabil- 
ity between the spin-down and spin-up states. Whereas, 
pure dephasing is preferred if the transition probability 
approaches zero. Hence, our strategy here is to relate 
the changes in the spin transition probability with the 
IET, and subsequently to the inelastic or pure dephasing 
mechanism. Apparently, we cannot use this strategy to 
calculate the dephasing rates. However, the changes in 
the transition probability can be related to the dephas- 
ing rates. This means that, if the transition probability 
approaches one, then the dominant dephasing rate is de- 
termined by the inelastic dephasing. 

In contrast, in the event where the transition probabil- 
ity approaches zero, then the pure dephasing will act as 
the dominant mechanism and determines the dephasing 
rate. Importantly, the ionization energy approximation 
requires that all the electrons are correlated and they 
are not free-electrons, be it in the conduction or valence 
band. Thus, this approach is not suitable for free-electron 
metallic substrates. Recall that our strategy here is to 
develop a theory to understand why the experimental 
data in Ref. Q is different from Ref. 



A. Dephasing rates 

First, the total dephasing rate (relaxation lifetime) for 
QDs can be written as [HJ 



1 



1 



rQD 



' imp 



r c-c T c-ph(SUB) 

1 1 

+- + , 

7"def Tspin 



~e-ph(QD) 



(22) 



where, r^ D is the total dephasing rate in a QD, while 
the other respective dephasing rates are due to e(QD)- 
e(QD), e(QD)-p/i(SUB) and e(QD)-p/i(QD) interactions, 
impurities, defects and spin. For non-magnetic and/or 
magnetic system above the Kondo temperature and in 
the absence of applied magnetic field, we can ignore 
Tspin, while Tj mp and rjof can be re-labeled as td that 
includes all types of defects (vacancy, interstitial, sub- 
stitution at different sites), including different types of 
elements. Thus, Eq. (|2"2"|) can now be rewritten as 
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T QD To-c T"c-ph(SUB) T"c-ph(QD) TD 



The next step is to invoke Eqs. (fl4|) and (|T5|) . which im- 
ply that the electrons in the QD (e(QD)) are influenced 
by the phonons from the QD (ph(QD)) and substrate 
(ph(SUB)). Therefore, 7 
can be simply written as 



b-c To-ph(suB) an d T _ ph (Q D ) 



T QD c(QD) 1 c(QD) 1 

T ph(SUB) T ph(QD) 



(24) 



The reasoning used to arrive at Eq. implies that 
the effects of defects, impurities, and e-e interaction in 
the QDs and substrate have been captured via e(QD), 
ph(SUB) and ph(QD). In other words, we can write 
Eq. (|2"2"|) in the form of Eq. (|2"4"]) because any systematic 
changes to these effects also affect the valence states (of 
each element) and/or the ionization energy of the system 
systematically. All our discussion in the previous sec- 
tions are based on Eq. (|24[) . For example, longer r^ is 
only possible if t SUB ^ D > rQ D (^ D > hw S p ^ QD ): adia- 
batically decoupled and pure dephasing as the dominant 
mechanism, or £SUB,qd < t qd ( ^qd < ^sub,q D) . adia _ 

batically coupled and inelastic dephasing as the dominant 
mechanism. 



B. Spin transition probability and the relaxation 
lifetime 

Claim 3: Even though we did not calculate the 
total dephasing rate in QDs explicitly under different 
perturbations (laser, electric and magnetic fields) using 
Eq. ([24]) . but we have invoked that the transition proba- 
bility of a confined electron in a QD from its ground to 
the excited state is proportional to the inverse external 
timescales, l/i SUB and l/i QD . For example, t SUB > r QD 
and tQ B > tQ d give rise to pure dephasing, while 
7j SUB < t qd and t QD < t qd lead to inelastic dephas- 
ing as discussed earlier. 

Proof 3: Here we will prove why such proportionality 
(transition probability oc l/t SUB >Q D ) is valid by means 
of the well known hamiltonian of an electron that starts 
out as a spin-up electron in the presence of a rotating (lo) 
magnetic field (-Bo) at an angle, a. The hamiltonian and 
its transition probability to spin down is given by [45[ 



W)=*£[A], 



(25) 



.4 



sin a cos(tut)(r x + sin a sm(ujt)<7y + cos aa z 



\(x(t)\x-(t))\ 2 = 







— sin a sin 









(26) 



A 



where, a x , a y and a z are the Pauli spin matrices, 
\Aj 2 + lo\ — 2ujuj\ cos a, and X±(A) denote the nor- 



— . (23) malized eigenspinors. Here, 



x(t) 



xe 



/„x (wi — W COS Oi) . 
cos(Ai/2) - '- sin(At/2) 



'*+(*)+* 



A 



sinasin(Ai/2) 



e^ X -(*),(27) 



X+(t) = (cos(a/2),e^sin(a/2)), 
X-{t) = {e~ wt sin(a/2),-cos(a/2)). 

The angular velocity, to = 1/T ex t refers to the charac- 
teristic time for the change in the Hamiltonian (external) 
given in Eq. (|25|) . while oj\ = eB /m = 1/T int refers to 
the characteristic time for the changes in the wave func- 
tion (internal). 



•i-Kext) > T(int): adiabatic 
— T(exl) < T(inl): non-adiabatic 

[(7Wr 6!tt ) sina] 2 [(T iat /T cxt )sin a\ 2 




2-nT 4ttT QttT KttT IOjtT t,me 

FIG. 4: Spin transition probability from spin-up to spin- 
down based on Eq. (|26|l . The label, T cx t > Tint indicates the 
system approximately satisfies the adiabatic condition, while 
Text < Tint implies the spins are non-adiabatically coupled 
to the external magnetic field (-Bo)- Thus, 1/T cxt determines 
the spin transition probability (spin-up to spin-down) . These 
plots have peaks at [(Text/Text) sina] 2 and [(Tint/Text) sina] 2 
for T ext < Tl n t and T cxt > Ti nt , respectively. The timescale, T 
labeled with 2nnT, n = 0, 1, in the time— axis corresponds 
to T ext and T int for T cxt < T int and T cxt > T int , respectively. 



The transition probability given in Eq. (|26|) is plotted 
in Fig.[¥]that captures the anticipated proportionality be- 
tween the transition probability and 1/T cxt . Recall here 
that the spin's timescale is T int . The plot labeled with 
Text > Ti n t in Fig. [4] implies adiabatic condition, in which, 
the spin transition probability is approximately indepen- 
dent of the external magnetic field, Bq (T ext ). Whereas, 
for T GXt < Ti n t, there is a non-adiabatic coupling be- 
tween the external Bq and the spin transition probabil- 
ity in the system. This in turn implies that l/T 0Xt cx 
\{x(t)\x-(t))\ 2 , as claimed earlier ■ 

In other words, we have a clear indication that smaller 
T oxt (compared to Tj nt ) leads to larger l/T oxt , which 
gives rise to larger spin transition probability (from spin- 
up to spin-down that corresponds to inelastic dephas- 
ing). On the other hand, large T ext will lead us to 
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smaller spin transition probability (this will enhance 
the pure dephasing mechanism). Using Eq. we 
can calculate the peak value for the plot, T oxt < T n t 
(non-adiabatic), which is given by, [(T ox t/T ext ) sin a] 2 . 
The peak value for the T oxt > Tint (adiabatic) plot 
is [(Ti n t /Text) sin a] 2 . Therefore, the values for these 
peaks satisfy [(T ex t/T oxt ) sin a] 2 > [(T in t/T cxt ) sin a} 2 
(see Fig.0}. 

Note here that the frequency for both plots are kept 
constant for graphical convenience by using the same 
magnitude for T cxt and T nt where the magnitudes of 
T oxt and T nt are interchanged when the inequality in 
T ox t < T n t is changed from < to > . In summary, we have 
shown that the transition probability is proportional to 
the external timescale 1/T cxt , which has been correctly 
invoked to discuss the experimental data in the earlier 
sections. For example, 1/T cx t (related to rotating (oj) 
magnetic field) can be identified as the 1/< SUB or l/i^ D , 
which are related to the w^ B and , respectively in 
Eq. (fT3|) . On the other hand, 1/T; n t (wi) is related to the 

1/r QD (u ,QD)_ 

VI. CONCLUSIONS 

In this work, several electronic and phononic issues in 
quantum dots arc discussed theoretically (as listed be- 
low), which also agree with recent experimental results. 

• The Green's and Gauss theorems have been invoked 
to derive the required electronic isolation of the quantum 
dots from the substrate in order to justify the electronic 
confinement. 

• Three conditions, (i), (ii) and (iii) were introduced 



to explain the physical mechanism required to isolate 
the quantum dots electronically. These conditions were 
used to discuss the phononic isolation that could occur 
in quantum dots, for different doping elements in QDs 
and substrates. 

• The ionization energy and its approximation have 
been employed to quantitatively explain how one could 
isolate the quantum dots from their substrate and/or ma- 
trix by varying their elemental compositions. 

• When the quantum dots are adiabatically isolated 
from the substrate, the pure dephasing process turns to 
play the leading role in determining its relaxation life- 
time. On the other hand, if the phonons from the sub- 
strate are non-adiabatically coupled to the quantum dots, 
then the inelastic dephasing becomes the dominant mech- 
anism. 

• All the theoretical results presented here, which are 
based on the ionization energy approximation and quan- 
tum adiabatic theorem agree well with the recently re- 
ported experimental and theoretical results. 

Finally, the ionization energy method presented here 
can be used to evaluate the types of quantum dots that 
one can grow, whether they are isolated with respect to 
substrate electrons and phonons. 

Acknowledgments 

A.D.A. would like to thank the School of Physics, Uni- 
versity of Sydney for the USIRS award and Micheal De- 
lanty for pointing out Refs. [H, H3| on QD shape effect. 
A.D.A. also thanks Kithriammah Soosay for the sup- 
port. K.O. acknowledges the partial support from the 
Australian Research Council (ARC) and the CSIRO. 



[1] P. Borri, W. Langbein, U. Woggon, V. Stavarache, D. 

Reuter, A. D. Wieck, Phys. Rev. B 71, 115328 (2005). 
[2] S. Sanguinetti, E. Poliani, M. Bonfanti, M. Guzzi, E. 

Grilli, M. Gurioli, N. Koguchi, Phys. Rev. B 73, 125342 

(2006). 

[3] J. L. Skinner, Ann. Rev. Phys. Chem. 39, 463 (1988). 
[4] E. A. Muljarov and R. Zimmermann, Phys. Rev. Lett. 

93, 237401 (2004). 
[5] E. A. Zibik, T. Grange, B. A. Carpenter, R. Ferreira, G. 

Bastard, N. Q. Vinh, P. J. Phillips, M. J. Steer, M. Hop- 

kinson, J. W. Cockburn, M. S. Skolnick, L. R. Wilson, 

Phys. Rev. B. 77, 041307 (2008). 
[6] A. K. Sood, S. Dattagupta, Pramana J. Phys. 17, 315 

(1981). 

[7] T. E. Hodgson, L. Viola, I. DAmico, Phys. Rev. B 78, 
165311 (2008). 

[8] V. M. Fomin, V. N. Gladilin, J. T. Devreese, E. P. 
Pokatilov, S. N. Balaban, S. N. Klimin, Phys. Rev. B 
57, 2415 (1998). 

[9] E. P. Pokatilov, S. N. Klimin, V. M. Fomin, J. T. De- 
vreese, F. W. Wise, Phys. Rev. B 65, 075316 (2002). 
[10] B. Krummheuer, V. M. Axt, T. Kuhn, Phys. Rev. B 65, 
195313 (2002). 



[11] S. V. Goupalov, R. A. Suris, P. Lavallard, D. S. Citrin, 
IEEE J. Sel. Top. Quantum Electron. 8, 1009 (2002). 

[12] T. Grange, R. Ferreira, G. Bastard, Phys. Rev. B 76, 
241304 (2007). 

[13] P. G. Klemens, Phys. Rev. 148, 845 (1966). 

[14] S. Barman, G. P. Srivastava, Phys. Rev. B 69, 235208 
(2004). 

[15] G. P. Srivastava, The physics of phonons (Hilger, Bristol, 
1990). 

[16] X. Q. Li, H. Nakayama, Y. Arakawa, Phys. Rev. B 59, 
5069 (1999). 

[17] O. Verzelen, R. Ferreira, G. Bastard, Phys. Rev. B 62, 
4809 (2000). 

[18] L. Jacak, J. Krasnyj, D. Jacak, P. Machnikowski, Phys. 

Rev. B 65, 113305 (2002). 
[19] F. Vallee, F. Bogani, Phys. Rev. B 43, 12049 (1991). 
[20] M. J. Seong, O. I. Micic, A. J. Nozik, A. Mascarenhas, 

H. M. Cheong, Appl. Phys. Lett. 82, 185 (2003). 
[21] Y. N. Hwang, S. Shin, H. L. Park, S. H. Park, U. Kim, 

H. S. Jeong, E. J. Shin, D. Kim, Phys. Rev. B 54, 15120 

(1996). 

[22] C. Q. Sun, Prog. Solid State Chem. 35, 1 (2007). 

[23] Y. Masumoto, M. Ikezawa, B. R. Hyun, K. Takemoto, 



11 



M. Furuya, Phys. Status Solidi B 224, 613 (2001). 
[24] D. Mogilevtsev, A.P. Nisovtsev, S. Kilin, S.B. Cavalcanti, 

H.S. Brandi, L.E. Oliveira, Physica E 40, 1487 (2008). 
[25] E. A. Zibik, T. Grange, B. A. Carpenter, N. E. Porter, R. 

Ferreira, G. Bastard, D. Stehr, S. Winnerl, M. Helm, H. 

Y. Liu, M. S. Skolnick, L. R. Wilson, Nature Materials 

8, 803 (2009). 

[26] A. Chernikov, S. Horst, S. W. Koch, S. Chatterjee, W. 

W. Ruhle, J. Sweet, B. Richards, J. Hendrickson, G. 

Khitrova, H. M. Gibbs, D. Litvinov, D. Gerthsen, M. 

Wegener, Solid S tate Commun. 149, 148 5 (2009). 
[27] A. D. Arulsamy , |arXiv:physics/0702232] /9 (submitted); 

larXiv:0807.0745l (submitted) . 
[28] A. D. Arulsamy, X. Y. Cui, C. Stampfl, K. Ratnavelu, 

Phys. Status Solidi B 246, 1060 (2009). 
[29] A. D. Arulsamy, Physica C 356, 62 (2001); Phys. Lett. 

A 300, 691 (2002); Phys. Lett. A 334, 413 (2005). 
[30] A. D. Arulsamy, M. Fronzi, Physica E 41, 74 (2008). 
[31] A. D. Arulsamy, K. Ostrikov, Phys. Lett. A 373, 2267 

(2009). 

[32] A. D. Arulsamy, K. Ostrikov, J. Supercond. Nov. Magn. 
(2009) DOI 10.1007/sl0948-009-0499-0. 



[33] A. D. Arulsamy, A. E. Rider, Q. J. Cheng, S. Xu, K. 
Ostrikov, J. Appl. Phys. 105, 094314 (2009). 

[34] M. Xu, H. Zhao, K. Ostrikov, M. Y. Duan, L. X. Xu, J. 
Appl. Phys. 105, 043708 (2009). 

[35] W. R. Smythe, Static and dynamic electricity (McGraw- 
Hill Book Company Inc., 1939, New York) 

[36] Y. Kayanuma, Phys. Rev. B 44, 13085 (1991). 

[37] S. Le Goff, B. Stebe, Phys. Rev. B 47, 1383 (1993). 

[38] J. Kumar, S. Kapoor, S. K. Gupta, P. K. Sen, Phys. Rev. 
B 74, 115326 (2006). 

[39] D.R. Santos Jr., F. Qu, A.M. Alcalde, P.C. Morais, Phys- 
ica E 26, 331 (2005). 

[40] K. Ostrikov, Rev. Mod. Phys. 77, 489 (2005). 

[41] K. Ostrikov, Vacuum 83, 4 (2008). 

[42] R. W. G. Wyckoff, Crystal structures, Vol. I and II (John 

Wiley & Sons New York, 1963, 1964). 
[43] M. J. Winter (http:/ /www.webelemen ts.comfr . 
[44] B. Dora, M. Gulacsi, Phys. Rev. B 78, 165111 (2008). 
[45] D. J. Griffiths, Introduction to Quantum Mechanics, 

(Prentice-Hall, New Jersey, 2005). 



